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We improve on a method to compute the fermion contribution to the vacuum polarization 
energy of string-hke configurations in a non-Abehan gauge theory. We estabhsh tlic new 
method by numerically verifying the invariance under (a subset of) local gauge transforma- 
tions. This also provides further support for the use of spectral methods to compute vacuum 
polarization energies in general. We confirm that the vacuum energy in the MS renormal- 
' ization scheme is tiny as compared to the mass of the fluctuating fermion field. Numerical 

results for the physical on-shell scheme are also presented. 

I. INTRODUCTION 

00 ' 

CN I The electro-weak sector of the standard model suggests the existence of extended flux-tubes called 

Z- or cosmic strings, which may have profound consequences for cosmological questions [H-Q]- 
These configurations are, however, not protected by any topological argument and thus are clas- 
, sically unstable 0, @]. To investigate whether quantum effects provide dynamical stabilization, 

Qj ' it is very important to compute the vacuum polarization energy AE of the string configuration 

i-^i reliably. In ref. [3] we have recently provided a proof-of-principle computation of AE in an SU{2)l 

gauge theory in three spatial dimensions^. That approach had the drawback that it required the 
I introduction of an auxiliary field at spatial infinity to make various components of the calculation 

ly-^ ■ well-defined. In the present letter we demonstrate that the formulation simplifies considerably in 

Q> , a suitable set of gauges. In particular, the use of an auxiliary field at infinity is avoided altogether. 

I The string configuration is translationally invariant along its symmetry axis (which we choose to 

be z), i.e. it only depends on the distance p from the axis and the corresponding azimuthal angle if. 
Finiteness of the classical energy (per unit length) requires that the string configuration must be 
Q I pure gauge at spatial infinity, which turns out to have a non-trivial angular dependence due to 

the winding of the string. As a consequence, gauge variant functionals of the string fields, such as 
Feynman diagrams, are ill-defined. This is the major obstacle for a straightforward application of 
^ I the spectral methods [l^ to compute the vacuum polarization energy of a string. In ref. j3] this 

^ ' obstacle was circumvented by the introduction of a return string that unwound the fields at spatial 

infinity. In numerical calculations this has the disadvantage that spatial infinity can only be reached 
by extrapolation to very extended profiles of the return string. These wide extensions induce large 
impact parameters so that channels with very large angular momenta must be considered. Here we 
argue that there are particular gauges in which the computation of AE does not require any return 
string. The litmus test then is to establish the invariance of AE under changes of parameters that 
classify these gauges. We present numerical evidence to confirm that this is indeed the case. Our 
finding gives further support for the use of spectral methods in general, as it proves the equality of 
gauge variant and divergent Feynman diagrams, and (equally gauge variant and divergent) terms 
in the Born series. 

In more detail, the string configuration consists of SU{2) vector and Higgs fields, and 
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See also ref. [7[ for a discussion of earlier attempts to estimate AE for string configurations, both in Abelian 

and non-Abelian theories. 
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respectively: 
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Here, we have used temporal gauge = and the SU{2) isospin structure is written in explicit 
matrix notation. Moreover, v is the (classical) vacuum expectation value of the Higgs field and g 
is the gauge coupling constant introduced in the string parameterization for later convenience. 

The configuration ([T|) is commonly called a Z-string, because the corresponding component 
Z ~ tr J (Wt^) exhibits the spatial dependence of an Abelian string. The radial functions fcip) 
and fnip) approach unity at spatial infinity while they vanish for p = 0. They are the typical 
profiles of the Nielson-Olesen type of string [2ol |. The angle G [0, vr] is a free parameter that 
determines the relative weight of the gauge and Higgs profiles; it also measures the fractional flux 
carried by the Z-string. 

We are mainly interested in the contribution from the fermion fluctuations to the vacuum 
polarization energy because it dominates the boson contribution when the number N of other 
internal degrees of freedom {e.g. color) becomes large. Motivated by the standard model we consider 
a non-Abelian gauge theory in which the gauge field only couples to left-handed fermions. The 
fermion-string interaction is then given by the Lagrangian 

= *i7/. (d^ - igW^") Pl* + ^iT/.^^'Pif ^' - / ^ (^Pr + <^^Pl) * , (2) 

where = ^ (1 =t 75) are projection operators on right- and left-handed components, respec- 
tively. The strength of the Higgs-fermion interaction is parameterized by the Yukawa coupling 
constant /, so that the fermions acquire the mass m = vf via spontaneous symmetry breaking. 



II. CALCULATIONAL TECHNIQUES 

We extract the Dirac Hamiltonian from eq. ([2]) and perform a local gauge transformation H 
WHU where 



cos{n{p) 

U = —iPLTiex.p {in ■ T ^) + Pr with h= | — sin(n93) | . 





(3) 



Here ^ = ^(p) is an arbitrary radial function that defines a subset of gauge transformations. The 
transformed Dirac Hamiltonian becomes 
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Hmt = in 
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We have made explicit the dependence on the angles = — £, and ^ via the isospin matrix 

(— sin(x) — icos(x) e*"'^\ 
icos{x)e ^ sm(a;) ) 

while the explicit matrices in eqs. and ([5]) act in spinor space. 

The key idea is now to impose the boundary conditions ^(0) = and ^{po) = ^i. Together with 
the boundary conditions for the physical profiles fc and fu this defines a well-behaved scattering 
problem for which a scattering matrix and, more generally, a Jost function for momenta in the 
upper half complex plane can be straightforwardly computed. Furthermore, the Born series to 
these scattering data can be constructed by iterating H\at- In contrast to the exact Jost function, 
the individual terms in this series are gauge dependent, i.e. they vary with i{p). After collecting 
these ingredients we proceed as in ref. |7|: 

1. In each angular momentum channel we evaluate the Jost function for imaginary momenta 
from the Dirac Hamiltonian, eq. ([5]). To this end we continue the momentum /c, that is 
conjugate to the radial coordinate p, analytically by substituting A; — )■ ir with r being a 
real variable. From the Jost function we then subtract its first and second orders of the 
corresponding Born series. This difference is summed over angular momenta. The analytic 
continuation to imaginary momenta is important because it allows us to exchange sums over 



angular momenta with momentum integrals l4] and implicitly accounts for the bound state 
contribution to /\E. 

We introduce a fake boson field whose second Born order approximation of the Jost function 
has the same divergence structure as the combined third and fourth order Born terms for 
the fermion problem. Again this quantity is summed over angular momenta. The fake boson 
method is a computational trick to circumvent the very cumbersome evaluation of third and 
fourth order Born terms and Feynman diagrams. This simplification has been established 
for purely logarithmic divergent contributions. 

We integrate the difference of the two functions constructed above over imaginary momenta r, 
weighted by a kinematical factor characteristic for string-like configurations that are trans- 



lationally invariant along a fixed direction [21[ . The value of this integral is the phase shift 
contribution, AEs. 

4. We add back the first and second order Born contributions in form of renormalized Feynman 
diagrams of identical order in Hi^f We call this piece AE'fd and discuss the details of the 
necessary counterterms further below. 

5. Finally, we add back AEb which is the renormalized second order fake boson Feynman 
diagram that corresponds to the subtraction under 2. It should be emphasized that the 
renormalization of AEb is accomplished by the counterterms in the fermion sector. 

In total, the fermion contribution to the renormalized vacuum polarization energy per unit 
length of the string reads 



AE = AEs + AEpB + AEb 



(7) 
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III. NUMERICAL RESULTS 

In addition to the angle , the parameterization of the string background with the above motivated 
boundary conditions introduces three width parameters, wh, wg and w^^, 

fnip) = l- e~^ , = 1 -e"(^) and ^(/o) = 6 

This parameterization guarantees that the interaction Hamihonian, i^int is well defined at p — )• 
and no 1/p type singularity is encountered. Obviously, the litmus test for our calculation is that 
the final result for AE must not depend on the scale introduced in the gauge transformation 
profile. In our numerical studies we always assume the special case n = 1. 

Numerically the most cumbersome quantity is the phase shift contribution AEs- For small 
values of the scale parameters wh and wg, in particular, the calculation of AEs for a single 
background configuration takes several days of CPU time on a modern desktop computer. In 
the treatment of ref. 0] at least as much time is consumed for each set of variational parameters 
that characterizes the auxiliary return string.^ Technically, we compute the momentum integral 
in AEs with the methods described above only up to a numerical cut-off Tmax- For r > Tmaxj 
we approximate the integrand by an inverse power-law. The numerical cost of this approach is 
determined by the smallest width parameter in the problem: the smaller this width, the larger we 
have to take Tmax for the power-law approximation to be accurate. Since a larger value for Tmax 
also entails that more angular momentum channels must be summed, the numerics become quickly 
expensive for small widths. From various integration methods and treatments of the contributions 
from large angular and linear momenta, we estimate an overall numerical accuracy of 1-2 %, where 
the upper limit mainly applies to small widths. 



(8) 



A. Verification of the method 

Before turning to the full string problem we note that the fake boson simplification introduces 
additional parameters into the numerical calculation. We have numerically verified that these 
parameters have no effect on the final result. 

To verify the method, we will establish the invariance of the vacuum polarization energy within 
the subset of gauge transformations obtained by varying the width of the gauge transformation 
profile (,{p). It is sufficient to consider the MS renormalization scheme because any other scheme 
differs by finite counterterms that are manifestly gauge invariant functionals of the background 
field. We augment the MS scheme by the no-tadpole condition which adds the counterterm 

(9) 

such that the local first order Feynman diagram is exactly canceled. The corresponding result is 
shown in table U] for a typical set of parameters. Since we measure all energies in units of the 
fermion mass m = vf and all lengths in its inverse, in the MS scheme AE only depends on the 
specific shape of the background profiles. 

We observe that the variation of the total result, AE, with is significantly less than the 
estimated numerical error for AEs, even though some components of AE change by almost an 



^ The calculation in ref. 0] needs to be redone several times with varying sets of return string parameters in order 
to extrapolate to an infinitely distant return string. 



/:3 = ^tr, 
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A En 


AE 


0.5 


-0.2515 


0.3489 


0.0046 


0.1020 


1.0 


-0.0655 


0.1606 


0.0032 


0.0983 


2.0 


-0.0358 


0.1294 


0.0038 


0.0974 


3.0 


-0.0320 


0.1235 


0.0056 


0.0971 


4.0 


-0.0302 


0.1193 


0.0080 


0.0971 



TABLE I: Independence on the scale of the gauge transformation parameter. 
The other parameters are wh = wq = 2 and = 0.47r, i.e. the gauge field 
is fairly strong. All energies are given in units of the classical fermion mass 
m — vf. 



order of magnitude in the considered range'^ of w^. We find similar results for other variational 
parameters wh, wg and ^i. Within the numerical precision this confirms the gauge invariance 
of the vacuum polarization energy, at least for the subset of gauge transformation that we have 
tested. 

We have also verified that — )• vr— leaves AE unchanged within the numerical precision. This 
symmetry follows from the fact that for the fields, eqs. ([T|) this transformation equals a rotation by 
TT about the z-axis in isospace. However, acting with this rotation on the gauge transformation U 
in eq. ^ gives a completely different radial function ^(p) with appropriately modified boundary 
values. The exact Jost function remains unchanged while neither the Born terms nor the Feynman 
diagrams are separately invariant, only their combination is. 

It should be emphasized again that in the course of this computation, we have added and 
subtracted formally identical quantities that are per se divergent and gauge variant. Hence our 
study also confirms that their finite pieces are identical, an assertion that is vital for the use of 
our spectral methods. Previously, such an identity had only be shown for the leading order of the 



Born and Feynman series within dimensional regularization [19|]. 

The numerical data also confirms our previous result [3] that AE is very small (as compared 
e.g. against the the fermion mass m) within the MS scheme. Our previous findings were, however, 
less accurate since they also required an extrapolation to an infinitely distant return string. 

In the left panel of figure [1] we display the dependence of AE (in the MS scheme) on the angle 
^1 that characterizes the relative strength of the gauge field and Higgs background. While A£'(^i) 
is monotonously increasing with (i.e. with stronger gauge fields) it develops a minimum around 
^1 = 7r/4 when the width of the gauge field background is small. 



B. On— shell renormalization 



To discuss physical implications we need to impose the on-shell renormalization scheme. For an 
SU{2)l gauge theory the on-shell renormalization conditions and the corresponding determination 



of the counterterm coefficients have been discussed in ref. 2j].'^ To pass from MS to on-shell, we 



^ Considering even smaller values for becomes numerically even more expensive because the asymptotic behavior 

of the integrand for AEs sets in at momenta roughly proportional to 
* Though the renormalized Feynman diagram is properly displayed in ref. the formula for C4 is missing an 

overall factor 1/2. 
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FIG. 1: The parameter dependence of the vacuum polarization energy. Left paneh MS scheme, right paneh 
on-sheU scheme. The dots denote the computed data and hnes are sphne interpolations. 



have to add the finite and manifestly gauge invariant counterterms 



C2 



£et = citr, [W^.W^''] + ^ in \{{d^ - igW^) {8^^ - igW^') $ 
where W^i, = 9[^W^] — ig[W^, Wy\ is the field strength tensor. 



(10) 



The on-shell renormalization condition implies that the pole of the Higgs propagator remains 



at the tree level mass, rrih 



m 



(0) 



, with unit residue. This fixes the coefficients C2 and C4 and 



ensures the usual one-particle interpretation of the states created by the asymptotic Higgs field. 
Furthermore, we also demand that the residue of the gauge field propagator (in unitary gauge) 
is unity, so that asymptotic VK-fields create one-particle 1^-boson states. This condition deter- 
mines ci. The position of the pole in the gauge boson propagator is then a prediction, i.e. the 
physical (on-shell) VF-boson mass receives radiative corrections. In our conventions (with all ener- 
gies measured in units of m = fv) we have / = so that = 2\plrr?GF makes contact to the 
standard model parameters. Using / = 0.9 and g = 0.7 approximately reproduces the top-quark 
and VF-boson masses. Furthermore we use ^ih = mh/m = v/\/2 ~ 0.8. The corresponding results 
for A.E are shown in the right panel of figure [TJ The additional counterterm contribution in the 
on-shell scheme increases /\E slightly making its little binding effect from the MS scheme even 
smaller. 



C. Total Energy 



So far we have not considered the leading contribution to the energy per unit length of the string, 
i.e. the classical energy [3] 



cl 



27r / pdp I n? sin^ ,^1 



Ik 
P 



+ 



pp2 



(1 - fcf 



+ 



p 4/2 



(11) 



where all quantities under the integral are dimensionless. Assuming that there are A'^ internal 
degrees of freedom, e.g. N = 3 for color, the total energy 



E = Eri + NAE 



(12) 
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FIG. 2: Total energy, eq. ([T2|) for = 3, in the on-shell scheme. Left paneh g ~ f — 5.0, right paneh 
g = f = 10.0. Note the difference in the scales on the vertical axes. 



will always be larger than NAE at the present level of approximation since Ec\ is positive definite. 
In order to allow for quantum stabilization, the classical energy must be comparable to AE, i.e. tiny 
as well. This is not the case for standard model motivated parameters, which give Ed ~ 10 m^. 
Eq. (fTTI) shows that small a Ed requires large coupling constants g and /, or equivalently large 
masses of the fluctuating fermion. To demonstrate this behavior we consider Ed + 3AE for g = 
f = 5.0 and g = f = 10.0 in figure O For fermion masses of order 1.5 TeV we indeed observe a 
small binding as < for narrow fields. However, this may merely reflect the onset of the Landau 



ghost [23|, 124]. It thus seems unlikely that the Z-string can be stabilized by fluctuating fermions 



without adding fermion charge. 



IV. CONCLUSION 

We have presented a reliable computation of the vacuum polarization energy that originates from 
fermion fluctuations about a cosmic string. The present method is significantly more efficient 
than the only one available so far because it makes redundant the introduction of an auxiliary 
field near spatial infinity. We have resolved the obstacles that stem from the non-trivial structure 
of the individual string fields at spatial infinity by choosing a subset of gauges for which the 
scattering problem is well-behaved. We have verified the novel method by establishing invariance 
with respect to gauge transformations within this subset. This is far from trivial because in the 
process of computation formally identical but divergent gauge variant quantities are added and 
subtracted. As an important side-product we have generated further support for the approach to 



compute vacuum polarization contributions to observables by spectral methods |19l |. 

Our extensive numerical investigations indicate that the fermion contribution to the vacuum 
polarization energy produces some binding but it is far too small to overcome the large classical 
energy and fully stabilize cosmic strings, at least for parameters that are motived from the standard 
model. 

Another stabilization scenario has been suggested in the D = 2 + 1 model of ref. Due 
to symmetry restoration in the core of the string, the Higgs condensate vanishes locally and a 
significant number of bound states can be induced. Population of these bound states may generate 
a charged object that is energetically favored against an equal number of free fermions with mass 
m = vf. We stress that this binding energy is of the same order in the ^-expansion as the part of 
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the vacuum polarization energy that we have computed here. Hence the present calculation is a 
necessary ingredient in a future study of the quantum stabilization of charged cosmic strings. This 
study will be subject of a future paper that will also serve to provide the details of the present 
computation. 



Acknowledgments 



We thank N. Graham and O. Schroder for valuable discussion in early stages of this project. 



[1] A. Vilenkin, Phys. Kept. 121 (1985) 263. 

[2] R. H. Brandenberger, Pramana 51 (1998) 191. 

[3] A. Achucarro and C. J. A. Martins. [arXiv:081 1.12771 [astro-ph]. 

[4] E. J. Copcland and T. W. B. Kibble, arXiv:0911.1345 [hep-th]. 

[5] A. Achucarro and T. Vachaspati, Phys. Rcpt. 327 (2000) 347. 

[6] Y. Nambu, Nucl. Phys. B 130 (1977) 505. 

[7] H. Weigel, M. Quandt, N. Graham and O. Schroder, Nucl. Phys. B 831 (2010) 306 

[8] M. Bordag and K. Kirsten, Phys. Rev. D 60 (1999) 105019. 

[9] M. Bordag and I. Drozdov, Phys. Rev. D 68 (2003) 065026. 

[10] K. Langfeld, L. Moyaerts and H. Gies, Nucl. Phys. B 646 (2002) 158. 

[11] N. Graham, V. Khemani, M. Quandt, O. Schroder and H. Weigel, Nucl. Phys. B 707 (2005) 233. 

[12] A. Alonso Izquierdo, W. Garcia Fuertes, M. de la Torre Mayado and J. Mateos Guilarte, Phys. Rev. D 

70 (2004) 061702, Phys. Rev. D 71 (2005) 125010. 

[13] N. Graham, M. Quandt, O. Schroder and H. Weigel, Nucl. Phys. B 758 (2006) 112. 

[14] O. Schroder, N. Graham, M. Quandt and H. Weigel, J. Phys. A 41 (2008) 164049. 

[15] M. Groves and W. B. Perkins, Nucl. Phys. B 573 (2000) 449. 

[16] M. Bordag, Phys. Rev. D 67, 065001 (2003). 

[17] J. Baacke and N. Kevlishvih, Phys. Rev. D 78 (2008) 085008. 

[18] S. C. Davis, W. B. Perkins and A. C. Davis, Phys. Rev. D 62 (2000) 043503. 

[19] N. Graham, M. Quandt and H. Weigel, Lcct. Notes Phys. 777 (2009) 1. 

[20] H. B. Nielsen and P. Olesen, Nucl. Phys. B 61 (1973) 45. 

[21] N. Graham, R. L. Jaffe, M. Quandt and H. Weigel, Phys. Rev. Lett. 87 (2001) 131601. 

[22] E. Farhi, N. Graham, R. L. Jaffe, V. Khemani and H. Weigel, Nucl. Phys. B 665 (2003) 623. 

[23] G. Ripka and S. Kahana, Phys. Rev. D 36 (1987) 1233. 

[24] J. Hartmann, F. Beck and W. Bentz, Phys. Rev. C 50 (1994) 3088. 



